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Abstract 

We compute the couplings DD*tt,DiD*tt,D* D*ir, DiDitt using QCD sum rules. These couplings 
are important inputs in the meson exchange model calculations used to estimate the amount of 
J/ip absorption due to pions and rho mesons in heavy ion collisions. Our sum rules are constructed 



> 

O I at the first order in the pion momentum p^, which give the couplings that are not trivially related 



to the soft-pion theorem. Our calculated couplings, which somewhat depend upon the values of the 
heavy meson decay constants, are goD^-w = 8.2 it 0.1, goiD'-K = 15.8 ± 2 GeV, go'D'-K = 0.3 it 0.03 



> , 

, and 5ri:iiZ)i7r = 0.17 ± 0.04. 

00 
O 



PACS numbers: 11.55.Hx,13.20.Fc,12.38.Lg,14.40.Lb 



Electronic address: hung@phya. yonsei. ac.kr 



I. INTRODUCTION 



J /if) suppression [|| seems to be one of the most promising signal for QGP formation in 
RHIC. Indeed the recent data at CERN show an anomalous suppression of J / ip formation, 
which seems to be a consequence of QGP formation [Q]. However, before coming to a 
conclusion, one has to estimate the amount of J/ip suppression due to hadronic final state 
interactions. Consequently, there have been a number of works, using various models 0, ^ |^, 
1^, ^, ^, |TD|, |TT], |T^, |TB[, calculating the J/ip absorption cross section by light mesons. However, 
the estimate varies by an order of magnitude near threshold. At this stage, it is necessary 
to probe each model calculations further to spell out their corrections and uncertainties. 

In the effective meson-exchange model approaches [^, |10], |lT], 0, important ingredi- 
ents are the couplings between open charm mesons and light mesons. Precise determination 
of the couplings reduces uncertainties in the calculation of the dissociation process. More- 
over, a complete set of low-lying open charm mesons D(1870), i5*(2010), Di(2420) has to 
be included. This is especially necessary in order to probe the cross section above the 
threshold. To provide the basic building blocks for such a model calculation, we use QCD 
sum rules [0, ^ and calculate the couplings 00*71 ^D* 0*71 ^DiDin ^DiD*^ . These can be 
used to improve the existing effective model calculations, which can be applied to future 
calculation |T^ of the J / ip dissociation process. 

At present, there are two approaches in the literature to calculate the coupling in the QCD 
sum rule approach, light-cone QCD sum rules (LCQSR) [|17|, 0, |T9[ relying on the operator 
product expansion near the light cone, and the conventional QCD sum rules |^D|, 



based on the short-distance expansion. Predictions from LCQSR heavily depend on the 
twist-2 pion wave function at the middle point, whose value, however, has been at the core 



of debates for a long time |23[. Furthermore, the duality issue in constructing the continuum 
needs to be carefully considered [1^. Instead, the conventional QCD sum rules 12^, |2l|, |2^ 
do not suffer from such an uncertainty as the QCD parameters appearing in this approach 
are determined from the low-energy theorems such as PC AC and the soft-pion theorem. 



Though QCD duality again needs to be applied carefully in these sum rules ^], the 
uncertainties from the QCD side can be substantially reduced. 

In this work, we provide a systematic approach to calculate the couplings using the 
conventional QCD sum rules. Our sum rules will be constructed at the first order in the 
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pion momentum 0{pfj). We improve the previous sum rule calculations of the DD*tt and 
D*D*n couplings |2^. In particular, QCD duality will be correctly applied according 

to the recent suggestion in Refs. |2^, |2^. Also in constructing a sum rule for DD*^, we 
advocate the use of a different structure function. We then construct similar sum rules for 
the couplings, DiDitt and DiD*7t, which may be important for future calculation of the 
J/ip dissociation process accompanying Di meson [^. 



II. OPE FOR GENERAL CORRELATION FUNCTION 



In this section, we schematically describe a general procedure to perform the operator 
product expansion (OPE) of the general correlation function with a pion, 

n(p,g) = I J dW«-^(0|T{J(x)ric(x),c(0)r2n(0)}|7r(p)) . (1) 

Fi and r2 denote specific gamma matrices corresponding to the coupling of concern. In 
later sections, we will use this general prescription, by simple replacements, to calculate the 
OPE for the correlation function of concern. For instance, to calculate the DD*n coupling, 
we will have Fi = 7^ and F2 = ^75. 

To calculate the coupling that is not trivially related by chiral symmetry, we will consider 
the correlation function at the first order of the pion momentum in its expansion. The 
soft-pion limit of the correlation function is just a chiral rotation of a vacuum correlation 
function {i.e., without the pion), which provides a coupling that is trivially related to the 
vacuum correlation function. For example, in the case of the pion-nucleon coupling ||2^, the 
soft-pion limit leads to Goldberger-Treiman relation with gj^ = 1, which gives the coupling 
30 % lower than its experimental value. To determine the coupling more precisely, one needs 
to go beyond the soft-pion limit. 

First, we restrict ourselves to the OPE diagram shown in Fig. |1| (a). For this diagram, 
we can rewrite the correlation function into the form 



(M 



rTr 



— 1^ + rric 



;T2Daa{k,p)Ti 



(2) 



(g — ky — ml 

The Roman indices denote colors. The free c-quark propagator has been used in obtaining 
this. Dab{k,p), which is shown by the blob in the figure, denotes the quark-antiquark 
component with a pion. This can be separated into three pieces depending on the Dirac 
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matrices involved, 



(3) 



The three invariant functions of k,p are defined by 

1 °° 1 ^ 
12 nl 

n=0 



idkof-^ • • • idka„ 

Qn 



12 ^ — ^ n\ 

n=0 

1 °° 1 ^ 

24 n! lokn, ■ ■ ■ lok, 

n=0 



idkdi^ ■ ■ ■ idka^ 

Qn 



in] 



5^'\k) , 
6^'\k) , 

d'-^^k) . (4) 



Since we are constructing the sum rules at the order 0{p^), we need to evaluate the pion 
matrix elements at 0{p^). Noting that the pion matrix elements involved are symmetric 
under exchanges of any pair of indices ai ■ ■ ■ a„ and using the soft-pion theorem and PCAC, 
we straightforwardly obtain the foUowings, 
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Here ttIq and 5^ are defined via 



Pa{9 9a2f3 + 9a3a29aif3) - (« ^ /3) .(5) 



{0\d9sG''^7f3uHp)) = tS'^Up'^ . 



(6) 



Up to twist-5, these are all the possibilities coming from the expansion of the quark-antiquark 
components at the order 0{p^). 

The additional contribution to the OPE is shown by Fig. |l] (b) where one gluon emitted 
from the c-quark propagator is combined with the quark-antiquark component. Specifically, 



the c-quark propagator with one gluon being attached is given by ||28 



2(A;2-m2)2 

where Qa/3 = t^Ga/s- The color matrices are normahzed via Tr(t^t^) = 6^^ /2. Taking the 
gluon stress tensor into the quark- ant iquark component, one can write down the correlation 
function into the form 

d'^k r 2(g - k)g'ys + - ^ + imc)i(ye5 



n(p,g)=2.y^Tr|. 



[{q — ky — m^]^ 



[j,a,xBP'^'\k,p) + 7-e^^"^D™^(A;,p)] Tij . (8) 
At the order the two functions appearing above are given by 



12 X 32 " " " ' ' dko 

Drap = -^^{Pagrp-ppgrx){27cY5^'\k) . (lO) 

Another function involving the pion matrix element of the form 

{Q\d'^->'^TgsQesu\'K{p)) (11) 

is of the order 0{mq), which therefore has been neglected. 

Once the Dirac matrices Fj (i = 1,2) are given, we can straightforwardly calculate the 
corresponding OPE from Eqs. (|^) and (^. These give the OPE up to twist-5 at the order 
0{pn). Below, we will use this formalism to calculate the couplings DD*7r, D*D*tt,DiDitt, 
and 0*0171 by choosing appropriate Dirac matrices for them. 

III. SUM RULE FOR DD*tt 

In this section, we construct a sum rule for the DD*n coupling using the correlation 
function 

U,{p,q) = t [ d'xe''^-^0\T{d{xMx),c{0)^l5umHp)) ■ (12) 



As two momenta are involved in the correlation function, one can separate the correlation 
function into the following two pieces 

^^. = Fi{p,q)p^ + F2{q,p)q^ . (13) 



We construct a sum rule for Fi{p = 0,q), which gives a couphng determined at the order 
0{pfj). Defining the D*Dti couphng by [[1] 



{D*{q)\D{q-p)Ti{p)) =gD*D^p-^ 



(14) 



and using 



(D|d75M|0) 
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{0\d'yf,c\D*) = mD*fD*e 



^l 1 



(15) 



the low-lying pole contribution to Fi{p = 0, q) is given by 

_ rn%mr,*fDfD*gDD*^ 



(16) 



A slightly different sum rule using the same correlation function can be found in Refs. ||20 
2l| . Specifically, the correlation function in that approach was decomposed into 



Ii^ = Ap^ + B{2q-p), 



(17) 



and a sum rule was constructed for the function A. Note, by comparing this with Eq. (P^, 
one can immediately see that A = F1 + F2/2. Thus, in the expansion in terms of the external 
momentum p^, the function A involves the term at the zeroth order in p^ [i.e. F2), which can 
be trivially obtained from a vacuum correlation function via the soft-pion theorem, as well 
as the term of the first order in p^ {i.e., Fi). To avoid the trivial contribution obtained from 
the soft-pion theorem, we choose to work with the decomposition of Eq. (0). Furthermore, 
as mentioned in Ref. |T^, the function A can contain some contribution from a possible 
resonance (scalar particle -Do), which can give an additional uncertainty in the prediction. 

Following the general strategy given in Sec. II, we obtain the OPE for the correlation 
function. 
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(18) 



6(g^ — m^)^ 3(g^ — m^)^ (g^ 

Note, the leading OPE has a simple pole in g^. According to QCD duality, higher resonance 
contributions lying along the positive are matched with the imaginary part of the OPE 
above a certain threshold Sq which is taken much higher than the low-lying pole. Since the 
simple pole structure l/(g^ — m^) does not have nonanalytic structures in the duality region 



{(f > 5'o) , it should not pick up the continuum contribution. However, it is an often practice 
that the double-variable dispersion relation is used to obtain a spectral density for a given 
OPE. Then by naively restricting the dispersion integral below the continuum threshold, 
one picks up the continuum contribution even from the OPE of the form l/{q^ — itl"^). 
In this prescription, the continuum contribution is a simple (and unphysical) pole at the 
continuum threshold [^, |2^. (See for example Eq.(3.12) of Ref. [^].) This pole at the 
continuum threshold does not resemble at all the higher resonance contributions lying along 
the positive g^. In fact, this pole at the continuum is mathematically spurious. To illustrate 
this in detail, let's determine the spectral density for the OPE l/(g^ — m^) from the double 
dispersion relation when the external momentum is zero. 

Under the successive Borel transformations ||2^, one can determine the spectral function 



b{s) = -e{s - mi) . (20) 

When we put it back to the double dispersion relation, we have to reproduce the OPE 
l/(g^ — m^). Anything additional to it is mathematically spurious. Using Eq.(PD|) in Eq.(|T^ 
and doing the integration by part, we obtain 
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(21) 



The second term is mathematically spurious as it is additional to the one that we had started 
with. But when it is restricted by the continuum threshold 5*0, the upper limit is changed 
to 5*0 and the second term has a pole at the continuum threshold. However, since the pole 
comes from the spurious term, its contribution to the sum rule is spurious. 

Nonetheless, one may argue from an intuition that the continuum contribution should be 
present as the current can couple to higher resonances. In fact, it may be possible to build 
such a contribution if one uses a more sophisticated current than the simple current of the 
form qVc. We believe that the absence of the continuum is due to a limitation of the current 
of the form qVc. But we believe that, as we demonstrated above, it is ad hoc to build the 
continuum contribution from the OPE of the form l/(g^ — m^). 



We now match Eq.(16) with Eq.(18) to get a sum rule for the D*Dtx coupling. For 



simplicity, we neglect the mass difference between D* and D and set them to mo* = tud 



rUav = {rriD* + mD)/2. Under the Borel transformation (with the Borel mass M^), the final 
sum rule reads, 
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(22) 



Here Ti denotes the transitions from the low-lying resonance to higher resonances. We will 
linearly fit the RHS within a Borel window to determine the coupling as well as the transition 
strength Ti. 



IV. SUM RULE FOR D*D*-k 



We now construct a sum rule for the 0*0*71 coupling. The D*D*tt sum rule can be 
constructed from the correlation function (by setting Fi = 7^ and r2 = 71/ in Eq. (|lD). 

^,.u{p,q) = t j rf'xe^''-"(0|T{(i(a;)7Mc(x),c(0)7.M(0)}|7r(p)) . (23) 

Saturating the correlation function by the D* intermediate state and introducing the cou- 
pling via 



{D*{q,e2)\'K{p)D*{q-p,ei)) = iY9D*D'ntaPt,ve1e^2P^q'' 
the low-lying pole contribution at the first order in is given by |]22| 
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The OPE part can be computed by following the general prescription described in Sec. II. 
After taking out the common factor of ^ap^yP^q'^ ■, we obtain the OPE side 
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(26) 



The terms involving rn^ are different from the one in Ref. p2[. By matching the two sides, 
we obtain the sum rule for the D*D*7t coupling 



,iml,-ml)/M^ 
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Here, T2 denotes the transitions from D* higher resonance states. 



(27) 



V. SUM RULE FOR DiDitt 

For the coupling DiDiir, we use the correlation function involving axial- vector currents, 
z [ d^xe^'^-^(0|T{d(x)7;.75c(x),c(0)7.75n(0)}|7r(p)) . (28) 



Comparing this correlation function with Eq. (P5D, one can easily see that the OPE for this 
correlation function can be obtained by replacing —ttLc in Eq. (^^. We introduce the 

DiDiTT coupling similarly as Eq. (P^. Then it is a simple matter to construct a sum rule 
for the DiDiTT coupling. Namely, by replacing rric —rric in Eq. (p7|), we have 



3 A 9' M2 9/^M4 

Again, T3 denotes the transitions from Di — > higher resonance states. 



(29) 



VI. SUM RULE FOR DiD*Tr 



The DiD*7T coupling can be calculated from the correlation function 

t j dW«-^'(0|T{J(x)7^c(x),c(0)7.75M(0)}|7r(p)) . (30) 

Refs. 0, |l^ calculated the coupling using the light-cone QCD sum rules. Both references 
considered the structure function proportional to g^y but the OPE in either approaches are 
different. Here we choose to work with the structure function proportional to g^pj^ + q^p^, 
which turns out to be rather simple. 



In constructing the phenomenological side, we follow Ref. where the structure pro- 
portional to g^pt/ + q^Pfj, is given by 

rriDi Id* foi 



9DiD*tt 



mD* q^ — mj-)* q^ — rn]-,^ 



(31) 



Note, gDiD*TT has one dimension due to the way that the coupling is introduced in Ref. |19 



which is in contrast to the other dimensionless couplings. The OPE side can be calculated 
straightforwardly. It takes the simple form of 



-A 



10 



ml 9 (g2 — m^)3 



(32) 
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The terms containing 5^ from the expansion of A{k,p) is canceled with the similar term 
coming from -D™/?- Setting ijid* = ttidi = ^ = {mD* +rnDi)/2 and matching the two sides, 
we obtain 



gDiD'Trfo'fDi + TiM"^ = f^e 



(m^~ml)/M'^ 



2- 

9M2 



(33) 



VII. ANALYSIS AND RESULTS 

In our analysis, we use the following set of the QCD parameters, 

ml = 0.8 GeV^ ; {qq) = (-0.24 GeV)^ ; S'^ = 0.2 GeV^ 



rrir 



1.34 GeV ; ^ = 131 MeV . 



For the hadron masses, we use 29 



rriD 



1.87 GeV ; mo* = 2.01 GeV ; mo, = 2.42 GeV . 



(34) 



(35) 



We plot the Borel curves for the couplings DD*^, DiD*^ in Figs. |^ and |^ using Eqs. (|22D 
and ( ^3]) respectively. The corresponding curves for the D*D*tx and DiDitt couplings are 
shown in Fig. |^. To get the couplings, we need to fit each curve with a straight line within 
an appropriately chosen Borel window. The intersection between the best fitting curve 
and the vertical axis at = gives the values fofD'gDD'n, fDjD*gDiD*7r, fh-gD'D'n 
and f£)-^gDiDiTT- All the curves are well-fitted with a straight line above minimum Borel 
mass, which depending on sum rules, ranging from 2-3 GeV^. In that region, the higher 
dimensional terms are suppressed. When we shift the Borel window by 0.5 GeV^, the results 
reduce by 10 %. 

Table | shows the best fit values and the chosen Borel window. For fDfD*gDD*-Ki our 
value is substantially smaller than 0.51 GeV^ obtained from the light-cone QCD sum rule 
analysis [0. The origin of the difference may be traced back to the use of the asymptotic 
pion wave functions. In the light-cone sum rules, the values of the twist-2 and twist-3 wave 
functions at the middle point, ^9,^(1/2) ~ 1.2 and ipp{l/2) ~ 1.5, enter as leading terms in 
the OPE. But in our sum rules, these middle points are replaced by the integrated strengths 
of the wave functions duipT,{u) = 1, du(pp{u) = 1, which are well- fixed by low-energy 
theorems. Though our approach is different from the similar calculations in Refs. ^0, 22 



our result for fnfD'gDD'n agrees with them. Also our value of f'^*gD*D*-K agrees with the 
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results in Ref. [|T^. Our result for fDifD*gDiD*n, 0.948 GeV^, is somewhat larger than 0.68 



Qgys ^Yom Ref.|T9|, which is obtained from the structure function proportional to g^^- 

To get the couplings, we need to determine f^, /d-^ and fa*. One may calculate these 
using QCD sum rules of two-point correlation function in the vacuum. Currently these 
values are not known precisely. According to Ref. they are fr, = nO MeV, fu* = 220 



MeV and Jdi = 240 MeV. A somewhat different set of the decay constants can be found in 
Ref.H, Id = 160 MeV, Jd* = 240 MeV and fo, = 300 MeV. Using these, we obtain the 
coupling constants, 

gDD*n = 8.29 (8.07) ; gD.n^n = 17.95 (13.6) GeV 

gD*D*. = 0.33 (0.278) ; gD,D,. = 0.208 (0.133) . (36) 

where the numbers (the ones in the parenthesis) are obtained by using the decay constants 
given in Ref. |TE[ (Ref. ||19||). To make a more precise prediction, one may need to determine 
the decay constants precisely. To summarize, we present in table || our results in comparison 
with the other previous calculations. 

These couplings are important ingredients for estimating the absorption cross section of 
J/ip hj IT mesons. Up to now meson exchange models in the calculation of the absorption 
cross section are based on a effective chiral lagrangian with only pseudoscalars {D, tt) and 
vector mesons {J/ip^D*). We believe that the addition of the axial partner of D* , namely 
the Di meson, will reduce the value for the existing calculation. The reason is the following. 
Consider the dissociation of J/ip by pions into D and D* . In the existing meson exchange 
calculations, the diagrams contributing to the process are the contributions from the s- 
channel D meson, the t-channel D* mesons and the direct four point coupling, which gives 
a non-trivial contribution. The form of the couplings are obtained from the chiral SU(4) 
lagrangian with vector mesons introduced in the massive Yang-Mills approach. We believe 
that the addition of Di meson will partly cancel the contribution from the direct four point 
coupling. This is so because this is precisely how the Adler consistency condition for the p—% 
forward scattering amplitude is obtained in the massive Yang Mills approach Namely, 
one has to introduce the axial partner of p, namely oi meson, whose contribution in the 
s-channel will cancel the direct four point coupling of the ppTrvr and make the amplitude 
vanish in the soft pion limit. Similar cancellation will occur between the direct four point 
coupling oiJ/ip — -K — D — D* and the J/ip + n^Di^D + D* contribution. Of course, 

11 



how big the cancellation actually is in this particular case has to be studied in detail and 
that will be done in the future work reported in Ref. |T6[ . 
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TABLE I: The best fitted values for the couphngs are listed here along with the chosen Borel 
window. Overal shift of the Borel window by 0.5 GeV^ to higher mass region gives 10 % error. 

Borel Window (GeV^) fitted value 

2.5 < M2 < 3.5 fofD'ODD*. = 0.31 GeV^ 

3 < M2 < 4 fDjD*9D,D*n = 0.948 GeV3 



2 < M2 < 3 fl*gD'D'n = 0.016 GeV^ 

2 < M2 < 3 /BifDiDiTT = 0.012 GeV^ 



TABLE II: Comparison of our results with the other published results. The results of Refs.[17, 18, 



19| are from light-cone QCD sum rules, the results from Refs.|20, p^] are from the conventional 



sum rules, and the result of Ref.|3C] is from the three-point sum rules 



9DiD*7T (GeV) 



9d*d*-k 



gOiDnr 



this work 



Ref.[|l7| 



Ref.lll 



Ref.[|T9| 
Refs.f^, |2 



Ref.[30| 



8.2 ±0.1 
12.5 ± 1 
11.85 ±2.1 

9±2 
5.7 ±4 



15.8 ±2 



10 ±2 



0.3 ± 0.03 



0.31 ±0.08 



0.35 ±0.08 



0.17 ±0.04 
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FIG. 1: The OPE diagrams considered in this work. The blob in (a) denotes the quark-antiquark 
component with a pion and the blob in (b) denotes quark-antiquark-gluon component with a pion. 
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FIG. 3: The Borel curves for the DiD*tt couphng from Eg. (|33|) . Here the RHS of the equation is 
plotted with respect to the Borel mass M^. 




M^(GeV^) 



FIG. 4: The Borel curves for the D*D*tt and DiDin couplings given in Eqs. (27) (p^). Here the 
RHS of the equations are plotted with respect to the Borel mass M^. 
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